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Abstract. Tunable interaction between two atoms in a cavity is realized by 
i-G ' interacting the two atoms with an extra controllable single-mode squeezed field. Such 

a controllable interaction can be further used to control entanglement between the 

two atoms against amplitude damping decoherence caused by spontaneous emissions. 

C^ ' For the independent amplitude damping decoherence channel, entanglement will be 

^ , , lost completely without controls, while it can be partially preserved by the proposed 

strategy. For the collective amplitude damping decoherence channel, our strategy can 
enhance the entanglement compared with the uncontrolled case when the entanglement 
of the uncontrolled stationary state is not too large. 
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1. Introduction 

Quantum entanglement [HElEllllElE] is a fundamental property of multi-body quantum 
systems that shows the non-local feature of quantum states. Quantum entanglement has 
been commonly recognized to be an essential physical resource in the implementation 
of high-speed quantum computation and high-security quantum communication. 

Many efforts have been made to create entanglement between decoupled quantum 
systems. One natural way is to introduce a simple intermediate device [3, El El [101 fTT] , 
e.g., a single- mode field or an additional particle, whose coherent interactions with the 
systems lead to their indirect interactions with each other. The intermediate device 
can also be measured to extract information about the quantum systems for quantum 
feedback controls [121 ESI Ej to manipulate the entanglement dynamics. One may 
also utilize a dissipative environment [151 tlHl [13 [E], e.g., a collective decoherence 
environment, to generate entanglement, interacted with which the system irreversibly 
decays to a stationary entangled state. 

However, in most circumstances, quantum entanglement tends to be destructed in 
environments [HI [20l [211 122]. For example, independent decoherence channels always 
lead to disentanglement [23] that is not recoverable by local operations and classical 
communications. 

Generally, non-local operations are required to effectively protect entanglement. 
However, a non-local Hamiltonian generated from the internal interaction between 
quantum systems, e.g., the dipole-dipole interaction between two atoms via the vacuum, 
is sometimes not a good choice, because disentanglement can also be induced by 
decoherence under these interactions. 

This paper introduces a single-mode squeezed field in a quantum cavity to realize 
non-local controllable interactions between two identical atoms in the weak coupling 
regime. By altering the parameter amplification coefficient of the squeezed field, one 
can continuously adjust the coupling strengths between atoms, which can be further used 
to control the final entanglement between the two atoms in presence of decoherence. It 
should be pointed out that there is another interesting work on coupling the two atoms 
via the squeezed vacuum [23]. Compared with the squeezed vacuum, the auxiliary 
squeezed field in the cavity is more controllable, which would be helpful to control the 
stationary concurrence. 

The paper is organized as follows: the physical model applied in the paper is 
formulated in Sec. [2l Entanglement control strategies are discussed for two-atom 
independent amplitude damping decoherence channels, collective amplitude damping 
decoherence channels, and their mixture, respectively in Sec. [3l [H and [5l Conclusions 
and a forecast of the future work are drawn in Sec. [61 
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2. Model Formulation 

Consider the system of two identical two-level atoms interacting with a squeezed single- 
mode field in a quantum cavity (see Fig. [1]). The total Hamiltonian of the atoms 
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Figure 1. Two atoms undergoing decoherence caused by spontaneous emissions interact with a 
single-mode squeezed field in a cavity. 



and the cavity mode can be described as below with h assumed to be 1 without loss of 
generality: 
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HAc = u^ca'a + ^Y.^f^ + Y.(' 



i^aa^ + e«*atai^): 



+ (ee-'^"a^' + eV"a'), (1) 

where the first two terms describe the free Hamiltonians of the cavity mode and the 
atoms; Uc is the frequency of the cavity mode and Ua is the inherent frequency of the 
atom corresponding to the energy separation between the ground state and the excited 
state of each atom; a is the annihilation operator of the cavity mode and az , i = 1,2 
is the 2;-axis Pauli operator of the i-th atom. The third term represents the interaction 
between the atoms and the cavity mode, in which a^ = ax ± icy , i = 1,2, are the 
ladder operators of the i-th atom. The complex coefficient 

S=fl.g{f^'y), 

is the inner product of the transition dipole moment fl of each atom and the coupling 
constant 
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where r^*-* is the position of the i-th atom; k and eg are the wave vector and unit 
polarization vector of the cavity mode; and V is the normalization volume of the 
cavity mode. The last term is the Hamiltonian of the squeezed cavity mode, where 
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the parameter amplification coefficient ^ and the frequency VL are continuously tunable. 
Such a manipulable standing squeezed field in a high-Q cavity is realizable by squeezed 
state engineering developed recently [261 I2Z|- Roughly speaking, a three-level atom 
in a ladder configuration is introduced to interact with the cavity mode. In addition, 
a classical field is used to manipulate the three-level atom, through which one can 
continuously adjust the squeezed coefficient ^ and the frequency Vt. 

In the weak coupling regime, i.e., A = cUq — tUc and |,^| ^ \^^^'^\i H^c can be 
diagonalized by the following unitary transform |28j : 



U = exp 



1 ^ 



e«aa« - e«*ata«) 



which, by taking the ffist-order approximation of e^^^A, gives the following expression: 
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where h.c. refers to Hermitian conjugate; 
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Further, by adiabatically eliminating the degrees of freedom of the cavity mode, 
the following reduced two-atom Hamiltonian can be obtained: 



Ha = yI1 ^'^ + (/^le-^^^*^^^ V«af + h.c.) + [^^,e~ 



-<^vi^)ai') + h.c. 



where the terms of individual atomic interaction with the cavity are omitted due to the 
fact that 

^»|e«|VA, |ee«|/A 

under the large detuning condition A ^ e^^\ Since the parameter amplification 
coefficient C, and the frequency Q are tunable parameters, we have two control parameters 
Hi and Q in H^- In the interaction picture, H^ can be expressed as: 



Hf 



/.ie-^^a«af +/.ie^'^^a«ai^M + Ue-^^^a«ai^) + /.^e^'^vi^Vf ) , (3) 



when the parameter Q is fixed to be 2uJa- 
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Besides the cavity mode, the atoms also interact with other modes in the 
environment, which leads to the atomic spontaneous emissions. In the case that the 
environmental modes are at the vacuum state, the dynamics of atoms can be described 
by the following master equation [HI [25] : 

p = -^[Hf + H,„ p]+J2 r., U^P^? - ^P^?V« - ^aJV^p) . (4) 
The parameters 



2^+ - 2 



are the spontaneous emission rates of the individual atoms, where p = \p\ is the 
magnitude of the transition dipole momentum, while 

ri2 = Tsi = rF( V12) (6) 

represent the collective spontaneous emission rates induced by the couphng between the 
atoms. The function F^kor^) can be expressed as [HI 



i^(Vi2) = ^<|(l-3cos2^ 



cos(A;o^i2) sin(A;on2) 
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+ sin^ 9 
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where fco = uja/c and 6 is the angle between the dipole moment vector p and the vector 
fi2 = f^^^ — r^^''; ri2 = |ri2| is the distance between the two atoms. The spontaneous 
emission process also introduces an additional coherent dipole-dipole interaction between 
the atoms: 

Tj I (1) (2) , (1) (2)x 

where the coefficient ri in Hyi can be written as [25J: 



r] = -T \[l — 3cos^ ' 



sin(A;o^i2) , cos(/i;ori2) 



^.^2gCOs(V,2), 

koru 



(Vl2)^ (koTu)^ 

3. Independent amplitude damping decoherence channel 

When the distance ri2 between the two atoms is far greater than the resonant wavelength 
l/fco = c/uja of the atom, i.e., fco'"i2 ^ 00, the amplitude damping decoherence of the 
two atoms can be taken independently. Consequently, from Eqs. (jH]) and ([7]), we have 
7], ri2, r2i — > 0, from which the following master equation holds: 

p = -i[Hf, p] + r©[ai^V + rD[ai'V, (§) 

where the superoperator D[L]p is defined as: 

D[L]p = LpLt - U^Lp - ^pLtL, 

and the two Lindblad terms D[crl ]p, D[o"_ ]p represent the amplitude damping 
decoherence channels acting on the two atoms with the damping rate F > 0. 
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To measure the quantum entanglement, we use the concurrence [3J between the two 
atoms of the quantum state p: 

C{p) = max{Ai - A2 - A3 - A4, 0}, (9) 

where \[s are the square roots of the eigenvalues, in decreasing order, of the matrix: 

and p* is the complex conjugate of p. 

It is known that, in absence of the squeezed field, a two- atom system will always 
be disentangled under independent amplitude damping decoherence channels (see, e.g., 
Ref. p3]), and this is not recoverable by any local operations. However, the entanglement 
can be partially protected via the intermediate squeezed field, because the solution p{t) 
of Eq. (IH]) tends to a stationary state 
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as a convex combination of a pure maximally entangled state 
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and a diagonal separable state 



Ps = diag(l-3/?, /5, /3, /5), 



where 



/3 



1 



1 



4pir 



The subscript "m" is an abbreviation of "maximally entangled" , and the subscript "s" 
refers to "separable" . The corresponding stationary concurrence is: 

2pi(r-pi) 



C{Po 



max 
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which, when the coupling strength pi is tuned to be 
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reaches its maximum value: 



Gn 



V5 + 1 
v/5-1 



0.31 >0. 



(12) 



The plot of C(poo) versus pi/F is shown in Fig. [2l 
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Figure 2. Plot of C(poo) versus fii/T. 



The cumbersome proof of Eqs. (TTOl) and (TT2|) is shown in Appendix A We adopt 
here an ideal model in which the two atoms have precise positions f^*-'. In real systems, 
position fluctuations are always presented, i.e.. 
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~f^' =f=<^) + 5r<*), i = l,2, 



~-li) . 



where r is the actual position of the i-th atom and Sf^"^^ is the corresponding 
fluctuation. From Eq. ([2]), the actual coupling coefficients should be 

^j(</.i+fc-(<5r<i)+5r<2))-) 
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^ifc-(<5r<i)+<5r<2)) ^i^i 
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which, consequently, fluctuates the phase by 

501 = 01 - 01 = fc . {6f^^^ + 6f^^^) = 6(j)^^^ + 6(f)^^^ (13) 

for the pure maximally entangled state pm (which now it should be Pm{<t'i) = Pm{4>i + 
5(pi)) due to the fluctuations of the positions of the atoms. Assume that S(f)^^^ and (50^^-* 
obey Gaussian distributions with means and variances 271 and 272, one can verify by 
averaging over the random fluctuations that the pure maximally entangled state pm is 
blurred into a mixed state: 
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Apparently, the resulting entanglement is also reduced. In fact, in this case, the 
stationary state should be 
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with a modified stationary concurrence 

The corresponding maximum stationary concurrence can be further calculated as: 

1 



Cn,ax = ^ (^V4e-(Ti+^2) + 1 _ 1 

Obviously, we have Cmax > 0, which means that our strategy is still valid compared 
with the case without the squeezed field. 

However, the maximum stationary concurrence is reduced by the dephasing effects 
caused by the fluctuations of the positions of the atoms. In order to estimate the 
influence of the fluctuations on the stationary entanglement, it can be estimated from 
Eq. ^ that: 

7j = -vari''*>^^*^' 



<i|^|\ar(|5r-<^)|)=27r2('- 



where A is the wavelength of the field in the cavity; var(5(/)*^*-') is the variance of (J^'-*-' 
and ^r*^*^ = a/vbj(|5^^*)|) represents the magnitude of the position fluctuation for the 
i-th atom. Therefore, if one is capable of trapping the atom in the cavity such that 

5r« < A, (14) 

the dephasing coefficients 7j can be neglected. This is possible under the present atom 
trapping and coohng technique since the wavelength A is of the order of /im (see, e.g., 
Refs. [291 SOI EH [32]). In this case, the perturbed maximum stationary concurrence 
Cmax is deviated slightly from the ideal maximum stationary concurrence Cmax (e.g., 
C'max/Cmax ~ 90% whcu Sr^'"'' / X = 0.05 as assumed in Ref. [29|). We can also see the 
influence from the following example with parameters given in Ref. [30J, in which the 
mass m of the atom (Cs atom), the oscillating frequency u of the external freedom of 
the atom (which is different from Ua), the effective temperature Tea of the atom, and 
the wavelength A of the field in the cavity are given as: 

m = 133 X 1.66 x lO^^^kg ^ 2.2 x lO'^^kg, 

w = 27r X 0.53 X lO^Hz ^ 3.3 x lO^Hz, 

Teff = 1.3 X 10-^K, A ^ 0.9 X lO^^m. 

Here, we choose an effective temperature Tes of the atom that is ten times greater than 
the lowest cooling temperature (13/iK) given in Ref. [30], under which the position of 
the atom can be taken as a classical parameter because 

1 
hhj < -ksTes, 

where ks is the Boltzmann constant. In this case, the position fluctuation of the atom 
can be estimated from 
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which leads to 

5r^^ = 
Thus, we have 



ksTcs 



1.38 X 10-23 X 1.3 X 10-4 



2.2 X 10-25 X 3.32 X 1012 
5r^^ 2.8 X 10-^ 



2.7 X lO^^m. 



0.03, 



A 0.9 X 10-6 
from which it can be calculated that C'max/C'max ~ 97%. 

Although the stationary concurrence may not be strong enough to be directly 
applied in quantum information processing and would be deteriorated by dephasing 
effects caused by noises such as the position fluctuations of the atoms, it is still hopeful 
to be used for entanglement protection. In fact, the fidelity between the stationary state 
Poo and the maximally entangled state pm can be calculated as: 

e-{7i+72)^^r - III 1 



F{poo) = tr(pooPn 

Optimally, it should be 
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(15) 



Since the maximum fidelity -Fmax is always larger than 0.5, we can, in principle, increase 
the stationary entanglement by introducing additional entanglement purification 
process [33l [3i] . 

To illustrate our proposal, let us discuss their applications in some typical 
circumstances. Firstly, consider the initial states pj at the maximally entangled state 
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and Po at the mixed entangled state 
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Moreover, let F = I/tq, where tq is the relaxing time constant. Simulation results are 
shown in Fig. [3l 

It is shown in Fig. [3] that the entanglement of the quantum states always decays 
to zero without control as what is known in the literature [23]. The corresponding 
stationary state is the two-atom ground state: 
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Figure 3. Plots of concurrence C{t) where the initial states are chosen as (a) maximally entangled 
state Pq and (b) mixed entangled state p^. The plus-sign lines denote the uncontrolled trajectories; 
the triangle lines are for the free trajectories in absence of control and decoherence; the solid lines 
represent the controlled trajectories with chosen parameters r = 2fj.iT/{4p1 + F^), which is the 
proportion of the maximally entangled state pm in Eq. pO|l . 



which is at the boundary of the set of all separable states. The superscript "m" refers 
to the "uncontrolled" system. When our strategy is applied, the entanglement can 
be remarkably retrieved against decoherence, and the maximum concurrence of the 
stationary state is 

v/5-1 



when 



L'mav — . ~ U.di, 



2/iir 1 
— — — = — ^ 45 
4/i? + r2 ^ ^-^^^ 



It is also noted that our strategy can enhance the entanglement of the stationary 
state of the naked atoms (i.e., neither control nor decoherence exist). 



4. Collective amplitude damping decoherence channel 

When the distance between the atoms is far shorter than the resonant wavelength of 
the atom, i.e., kori2 -^ 0, from Eqs. ([6]) and ([7]) we have: 
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which corresponds to a two-atom collective amplitude damping decoherence channel |35j . 
In this case, the master equation of the two atoms becomes: 
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Hf+r,^{a';'a'_>+a'ya 
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where the two-atom operator S'_ = cr_ + cr_ , and F > is the damping rate. Because 
the two atoms are very close to each other, from Eq. ([2]) the coupling strength between 
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each atom and the cavity can be taken as identical, i.e., e*^^^ = e*^^^ = e, so that the 
interaction Hamiltonian H'^ can be expressed as: 

Hf = ^^^ (e-*aW,f + e^^^ai^Vi^^) + /., (a^ai^^ + a«af ) , 

where 

/zie-^'^i = 2eeVA2, ^^, = [ep/A. 

In absence of the intermediate squeezed field, i.e., /ii = /i2 = 0, the stationary state 
of the two-atom system 

P^ = (1 - f^)pm + i^Po (17) 

is a convex combination of the maximally entangled state 
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and the two-atom ground state 
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where the weight n G [0, 1] is determined by the initial density matrix: 
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ft = tr 
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and 
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The resulting stationary concurrence is 

c(p:^) = 1 - ft. 

When the intermediate squeezed field is presented, the corresponding two-atom 
stationary state 

Poo = s~pm + rpm + (1 - s - r)ps (19) 

is a convex combination of the maximally entangled states pm and pm given in Eqs. ( TT8l) 
and (TTTj) respectively, and a diagonal separable state 

p, = diag iPi+P2,-- P2, - - /32, /32 - /3i j , 
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where 



2(r2 + 3/i?)' ^' 2(r2 + 3/^f) 

The weights s and r are, respectively, 

, 7 r^-Su? 2r/iiK 



6 r2 + 3/i2 ' r2 + 3/i2- 

It can be examined that when the parameter fii is in the range: 



V-9k2 + 22fi; - 12 /ii k + V-Qk^ + 22k - 12 



6 — Sk ~r~ 6 — Sk 

the resulting stationary concurrence C(poo) is superior to Ci^p^) without the 
intermediate squeezed field: 

(r2 + 2/i2 + 2r /ii)/s: 

r2 + 3/i2 

The interval given in (l20l) is nonempty if and only if 

11 1 



c(Poo)= " "„^:i7"' -i>c(p-j=i-K. (21) 



.13 < K < 1, (22) 

9 9' ^ ^ 

otherwise, our strategy is not capable of improving the stationary concurrence. 
Moreover, the maxima of C(poo) is achieved when 

2 



13 + 1 
and the corresponding maximum value is 



C -^^i^±^.-l 

'—'max „ rb ±. 

O 

The plots of the stationary concurrence versus the coupling strength for different k are 
shown in Fig. HI In Fig. H^, the controlled stationary concurrence is superior to the 
uncontrolled one only in the interval given by Eq. fl20l) . while, in Fig. |3)d, the controlled 
stationary concurrence is always better than the uncontrolled one. 

The fact that our strategy is effective only when the parameter k is sufficiently 
large comes from the competition between pm and pm in Eq. (IT9l) . where pm comes 
from the dissipation effect and pm is induced by our proposal. When n is close to 0, 
the dissipation dominates and hence the control fails, while, when n is close to 1, the 
control becomes effective. 

As has been indicated in Sec. |3l the fluctuations of the positions of the atoms 
would bring an uncertain phase shift for the maximally entangled state pm-, which may 
deteriorate the stationary concurrence. The calculations are like those in Sec. [3], so we 
omit them here. 

Fig. [5] shows some numerical examples, where the initial states are, respectively, 
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(a) (b) 




Figure 4. Plots of C(poo) versus /ii/F for (a) k — 0.85, (b) k = 1. The asterisk line is for the 
controlled stationary concurrence and the solid line is for the uncontrolled stationary concurrence. 
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Figure 5. Plots of the concurrence C(t) where the system state is initialized for (a) p\ and (b) p\. 
The plus-sign lines denote the uncontrolled trajectories; the triangle lines are for the free trajectories 
in absence of control and decoherence; the solid lines represent the controlled trajectories with 
different parameters r — 2r/iiK/(r^ + 3/if), which is the proportion of the squeezed field induced 
maximally entangled state p„i in Eq. (|19p . 



The simulation results show that the stationary state of the uncontrolled system 
under the collective amplitude damping decoherence may remain entangled which is 
quite different from the independent decoherence channel, and this feature has been 
utilized in the literature [UIIISIIIZ] to create entanglement between qubits. Our strategy 
may further increase the entanglement in the stationary entangled state, as shown in 
comparison between the plus-sign lines and the solid lines. 
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Another feature of the collective amplitude damping decoherence channel observed 
from Fig. [5] is that the maximum concurrence depends on the initial state. For certain 
values of r, our control strategy may have worse performance than that induced by 
the natural dissipation. Both Fig. [5^ and [5b provide such a case where the solid line 
(controlled trajectory) goes below the plus-sign line (uncontrolled trajectory). The 
corresponding control parameter /xi is outside the interval given in Eq. fl20|) . 

5. Mixed amplitude damping decoherence channel 

In actual experiments, the decoherence channel is never perfectly collective, because it 
is hard to place the two atoms in an cavity close enough. The existing atom trapping 
and cooling techniques [36l |37j can only hold two atoms approximately at the distance 
of the same order of the resonant wavelength of the atom. Thus, it is more realistic 
to treat the resulting decoherence channel as a mixture of an independent amplitude 
damping decoherence channel and a collective amplitude damping decoherence channel, 
as shown in the following master equation: 

p = - ^[Hf,p] + J:rV[a^^]p + F,, (a«paf - i{af a«,p} 

+ r,JaiV?)-l{aiV^p}V (23) 



2 

where < F12 < F. 

It can be verified that the stationary state of the uncontrolled system is nothing but 
the separable two-atom ground state p^ = |00)(00| in which entanglement completely 
disappears, as well as in the case of the independent decoherence channel. By 
introducing the intermediate squeezed field, we can stabilize the system at the same 
stationary state given in Eq. (llOj) . 

6. Conclusions 

In summary, we proposed a two-atom entanglement control strategy, via a controllable 
squeezed field coupled to the two atoms, to protect entanglement from the spontaneous 
emission process. The parameter amplification coefficient of the squeezed field can be 
tuned to generate a non-local Hamiltonian, which can be used to maintain entanglement 
of the two-atom states against decoherence. For the independent amplitude damping 
decoherence channel, we can partially recover the entanglement of the quantum 
state which otherwise will be completely lost. For the collective amplitude damping 
decoherence channel, our strategy can effectively enhance the entanglement of the 
stationary state compared with the dissipation-induced strategies provided that the 
uncontrolled stationary state is not tightly entangled. 

The proposed entanglement control strategy is an open-loop control strategy, where 
no measurements are done during the course of control. Such control strategies require 
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exact values of the system parameters, and can badly suffer from the uncertainty of 
these parameters, which may bring remarkable derivation of the stationary concurrence 
from the ideal values. This problem is hopefully solvable by quantum feedback controls. 
Another direction of the succeeding research will be the application in solid 
state systems. In such systems, controllable coupling between qubits is easier to be 
achieved |3Hl [39] compared with the optical systems. However, interactions between 
the solid-state systems and their environments are more complicated, which may lead 
to non-Markovian noises |10]. To what extent the controllable non-local unitary opera- 
tions can preserve entanglement against non-Markovian noises is an interesting problem 
to be explored, for which existing decoherence suppression strategies [HI [121 [33] may 
be helpful. 



Appendix A. Proof of the Equations ( ITOl l and ( [T2ll 



Firstly, we transform the control model ([8]) from the complex matrix space into the real 
vector space, i.e., the so-called coherence vector picture [HI [13 [16], [U]. With respect 
to the inner product {X,Y) = tr(X"''y), we define the following orthonormal basis for 
all two-atom operators: 



where 



f-L . O^ ny Q^ ny 1^(1) -a^^^ -a^^'^ i^x^^) -a^^)a 

S 2-'4x4? "14) "145 "23) "23) 9 a; ) 2 V ' O^ ' O « ' 9 
1 (1) (2) 1„(1)^(2) 1^(1)^(2) Q^ Q. 1^(1)^(2)1 

2*^2 "^x ) 2 y 2 ' 2 y ' i4)"23)2^2 ^z f ) 
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z ) 



/ 



O^ — 



QX 

^^23 



V2 



l\ 



o\ 



/ 



Q.y — -— 



( 



) '^'23 



v^ 



0\ 



(A.l) 



^ 



14 



V^ 



/l 



V 







\ 



-1 



) '^'23 



V2 



/o 



V 



Under this basis, the system density matrix can be expressed as: 

15 



P — T-^4x4 



y^r»i^i 



i=l 



\ 
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where rui = tT{Qip) and f2j, i = 1, ■ ■ ■ , 15, are the basis matrices in Eq. (lA.ip except 
^/4x4- rn = (mi, ■ ■ ■ , mis)-^ is called the coherence vector of p. 

In the coherence vector picture, the master equation (JHl) can be rewritten 
as miSlllSlllTl: 

m = OAfTi + Dm + g, (A. 2) 

where the orthogonal matrix Oa is the adjoint representation [IS] of —iH^. The afiine 
term "Dm + g" is that of the Lindblad terms: 



.(1) 



rD[(Tr> + r©[(Tr> 



(2)i 



where D < and g is a constant vector. Further, divide m into the following sub- 
vectors: 



y \T 



where 



mP = (mi4,m5'4,m23,m|3) 

m'^ = (mi4,m23,m2^)^, 

m^ = {mM,myo,mox,moy,m^z,m^^,my^,m^yf, 



mi4 = tr(n°4p), m^3 = tr^Q^^p), a,p = x, y, z, 



(A.3) 



TUij = tr 



^ (1) (2) , 

2^i ^} ) P 



hj = 0,x,y,z, 



and aQ =/2x2, "^ = 1,2 are the 2x2 identity operators acting on the i-th. atom. Then, 
(]A.2|) can be grouped into: 

4 

mP = ^UiOlm'' \ DPm\ 



where 



4 



m'' 



■nf 



^ M,(-0,''^)mP + D'^m'' + g^, 

4 



(A.4) 
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Ml = OyUi COS 01, U2 = O/ii Sm(/)i, Us = 0/i2 COS 02, M4 
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1 \ / 2V2T \ 



D"^ = ~AT \ 10 

-V2 2 



, 9- 



V 








(A.5) 



/ 



Of are all skew-symmetric matrices and D"^ < 0, thereby 

— {{m'fm') = {m'fD'm' < 0, Vm' ^ 0, 

which implies that m'^ ^ when t — >■ 00. 

With simple calculations, the following stationary solution can be obtained for Eq. 

rrf = 0, m23(oo) = m%{(X)) = m^^ioo) = 0, 
V2fiiT 



m^^[oo) 



ml^{oo) 



m.Aoo 



TC 



^/ilF . / TV 



V2: 



m.Aoo] 



^2^ 



from which we can obtain the corresponding decomposition (fTOj) of the stationary state 

Poo- 

Going back to the density matrix, one can find that the stationary state has the 
following form: 

I a w \ 

b z 
z* c 
w* d 

whose concurrence can be analytically solved to be [191 EDI l2T] : 

C{poo) = 2max{|w| — v6c, \z\ — vofi, 0}. (A. 7) 

The above equation leads to the stationary concurrence C(poo) given in Eq. ( fT2l) . 



Pc 



v 



(A.6) 
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Appendix B. Proof of the Equations (TT71) - (pT1 ) 



Similar to what we have done in Appendix A, the controlled master equation (1161) can 
be grouped into the following control equations in the coherence vector picture: 

8/ii sin Oiml^ — AVml^, 

8/ii cos6'imf4 - 8/ii sin6'imi4 - 4rmi4 + 2^/21 + 4Tm2^, 
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(B.l) 



Using a squeezed field to protect two-atom entanglement against spontaneous emissionslS 

As well as the independent amplitude damping decoherence model, the sub-vector 
m*^ = [m^Q, ■ ■ ■ , rn^y)'^ always goes to zero when t —>■ oo, and it will not affect C{poq), so 
we will not discuss m'^{t) here. The fourth and the last equations in Eq. ( JB.ll) implies a 
conservation law: 

3 

m,,{t) + V2ml^{t) = m^^(to) + V2m^^{to) = 2k - -. 

Substituting /ii = /i2 = into Eq. ( IB. II) . we have the following uncontrolled stationary 
solution: 

m^4(cx3) = m^l^^i^oo) = ml^{oo) = m23(oo) = 0, 

m^4(oo) = — K, m23(oo) = — (k- 1), m^^(oo) = k - -, 

from which the decomposition flTTl) can be obtained and C{p'^) = 1 — n. 
Further, we can obtain the controlled stationary solution of Eq. f IB.ip : 

m,4(oo) = ^^-^^^ cos (01 - - j , m,,{oo) = ^;^^^ sm (0, - - 



m|3(oo) = m^3(oo) = 0, ml^{oo) = ^.^ ^ ^y 

1 pIk X I \ ^ ( 1 , /^i'^ 

which leads to the decomposition flT9|) . Further, it can be calculated that 



C(poo) = max{Fi(/ii),F2(/i2),0}, 

where 

F,(^,) = -. + 1 + ^^^-:^n - P + 3;.^ "' 

(r^ + 2/if + 2r/ii) 

i^2(/.i) = ^^^^^^^ «:-!. 

It is easy to verify that -Fi(/xi) monotonically decreases when the control parameter \i\ 
increases, thereby 

Also, we can obtain that 

F2(/ii)>l-«: = C(p^)>0, 
when \ix satisfies Eq. fl20l) . In conclusion, we arrive at 
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